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Abstract. The temperature of hot accretion flows around black holes is sufficiently high for the ignition of 
nuclear reactions. This is potentially an important nucleosynthesis mechanism in the universe. As the first step in 
studying this problem, we need to measure physical quantities such as density and temperature of the accretion 
flow. In usual studies of the hot accretion flow, viscous dissipation is considered to be the only heating mechanism, 
while the heating caused by nuclear reactions is not considered. In this paper, we investigate whether the energy 
generation rate of nuclear reaction is important compared to the viscous heating. Our calculation indicates that 
the former is at most one percent of the latter and thus is not important. The dynamics of accretion flow can be 
therefore calculated in the usual way, without the need to consider heating due to nuclear reactions. 
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1. Introduction 

Observations imply that there maybe a supermassive 
black hole at the center of almost every galaxies. There 
is a large fraction for which the accretion flow can be 
described by a hot accretion flow, such as advection- 
dominated accretion flow (ADAF; Narayan & Yi 1994, 
1995; Abramowicz et al. 1995; for a review of ADAF, see 
Narayan, Mahadevan & Quataert 1998; for a review of the 
applications of ADAF, see Yuan 2007) and luminous hot 
accretion flow (LHAF; Yuan 2001). Hot accretion flow also 
exists in some states of black hole X-ray binaries. The tem- 
perature of the hot flow is very high, TJ ~ 10 12 K/r, r being 
the radius in unit of Schwarzschild radius r s . This temper- 
ature is high enough to ignite significant nuclear reactions 
in the hot flow. This process is potentially interesting be- 
cause it is a promising means of producing new elements. 
The analysis of observations by Juarez et al. (2009; see 
also Hamann & Ferland 1999) found that the metallicity 
of the BLR of redshift quasars is very high. One expla- 
nation is that the massive stars had formed before the 
quasars were created, and these stars then enriched the 
metallicities of quasars. It is also possible that BLR gas 
originates in accretion flows. Although accretion flows in 
quasars are generally assumed to be in the regime of higher 



accretion rates than that of a hot accretion flow, i.e., the 
standard thin disk, the detailed structure of the accretion 
flow is still largely unknown because the thin disk can not 
explain the origin of the X-ray emission widely observed 
in quasars. The presence of X-ray emission implies that 
ADAF or LHAF-like hot accretion flow is also likely to 
exist in quasars, not only in some low-luminosity AGNs. 
One possible reason for the high metallicity in quasars is 
therefore that nucleosynthesis occurs in the accretion flow 
(see the review paper of Hamann & Ferland 1999). 

Some people have already studied the nuclear reactions 
and nucleosynthesis in accretion flows around black holes 
(e.g., Jin et al. 1989, Arai & Hashimoto 1992, Chakrabarti 
& Mukhopadhyay 1999, Mukhopadhyay & Chakrabarti 
2000, Mukhopadhyay & Chakrabarti 2001, Hu & Peng 
2008). In these works, however, the accretion flow models 
adopted are quite different from ADAFs. The works of Jin 
et al. and Arai & Hashimoto are based on a "thick disk" , 
while the two works by Mukhopadhyay & Chakrabarti are 
based on a simplified inviscid accretion flow, which dif- 
fers again from ADAFs. These particular accretion mod- 
els have received little attention these years. Hu & Peng 
(2008) studied nucleosynthesis in ADAFs, but their work 
was based on the self-similar solution of ADAF. This so- 
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lution is affected by significant error in the inner region of 
the accretion flow, where most nuclear reactions occur. 

We therefore revisit the study of nucleosynthesis in 
hot accretion flows. We first constrain the dynamics of 
accretion flows. Although the global solution of ADAF 
has been known for many years, since nuclear reactions 
have not previously been taken into account, we need to 
investigate whether they are important compared to the 
viscous heating. If they are then, the self-consistent global 
solution of ADAF needs to be recalculated. We neglect the 
photodisintcgration because the optical depth in ADAFs 
is very small and most of the photos can escape from the 
hot accretion flow without being scattered. 

The structure of the paper is as follows. We briefly in- 
troduce ADAFs in Sect. 2.1 and the calculation method of 
energy generation rate in Sect. 2.2. The calculation results 
will be shown in Sect. 2.3. In Sect. 3 we summarize our 
results. 

2. Energy generation rate in hot accretion flow 

2.1. Advection dominated accretion flow 

The dynamics of the ADAF is described by the continuity 
equation, the radial and azimuthal components of momen- 
tum equation, and the energy equations of electrons and 
ions (Narayan, Mahadevan & Quataert 1998) 
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where M out is the mass-flow rate at the outer boundary 
(-Rout) of the flow, the exponent s describes the strength 
of the outflow, S is the fraction of turbulent viscous en- 
ergy that directly heats the electrons, q + is the viscosity 
heating rate, qi e is the Coulomb energy exchange rate be- 
tween electrons and ions, and q~ is the radiative cooling 
rate per unit volume. The main radiative processes are 
synchrotron emission, bremsstrahlung, and their respec- 
tive Comptonizations. The eigenvalue j corresponds to the 
specific angular momentum loss at the inner boundary, £j 
(e e ) is the internal energy of ions (electrons), and Pi (P e ) 
is the pressure due to ions (electrons). 

The global solution of these equations needs to sat- 
isfy three conditions. The first one is the no-torque con- 
dition at the horizon. This condition can be automati- 
cally satisfied if we adopt the viscous description of T rr j, = 



ce{P gas + Pmag) , where i>0 is the viscous stress tensor, and 
Pgas and Pmag are the gas pressure and magnetic pressure, 
respectively. The second condition is the transonic condi- 
tion. When the matter enters into the horizon of the black 
hole, its speed should be equal to the light speed. A sonic 
point must therefore exist where the radial speed of flow 
equals the local sonic speed. The solution should finally 
satisfy the boundary condition at R ou t ■ Here we chose the 
temperature of ions and electrons T^ e , and the ratio of the 
radial velocity of the flows to the local speed of sound v/c s 
at Rout- We adopt the Paczynski & Wiita (1980) poten- 
tial to mimic the geometry of a Schwarzschild black hole. 
We adjust the eigenvalue j and use the shooting method 
to derive the global solution that satisfies the conditions 
above. For the parameters of the hot accretion flow, we 
choose their typical values constrained by modeling black 
hole X-ray binaries (e.g., Yuan, Cui & Narayan 2005) and 
the supermassive black hole in our Galactic center, Sgr 
A*, which is the most well studied supermassive black 
hole system observationally (Yuan, Quataert & Narayan 
2003). They are the viscous parameter a = 0.3 (but see 
discussions below), (i = 0.9, s = 0.3, and 5 = 0.5, where (3 
is the ratio of the gas pressure to the total pressure, s de- 
scribes the strength of the outflow, and 5 is the fraction of 
viscous dissipation that directly heats the electrons. The 
detailed numerical approach can be found in Yuan (200 f) 
or Yuan et al. (2003). 

As is well known, when the accretion rate is very low, 
almost all of the viscously dissipated energy is stored in 
the accretion flow and advected into the horizon of the 
black hole, rather than radiated away. This is because 
the radiative timescale is far longer than the accretion 
timcscale. However, when the accretion rate increases, 
the radiative timescale becomes shorter, and less energy 
is advected into the black hole. In the regime of LHAF, 
the energy advection even becomes negative. In this case, 
the radiative efficiency of the hot accretion flow increases, 
even reaching that of the standard thin disk (Shakura & 
Sunyaev 1973). Unlike the standard thin disk, ADAF is 
of course, still hot and geometrically thick. The tempera- 
ture of the ions is almost always virial, and is a power-law 
function of radius. 



2.2. Calculation of nuclear reaction heating 

We make the following assumptions. The first is that the 
element abundance in the accretion flow is the same as 
that in the solar atmosphere (Grevesse & Sauval 2000). 
The companion of the stellar mass black hole could be a 
main-sequence star, and the accretion matter of a super- 
massive black hole at the galactic center is presumably 
supplied by a number of nearby stars. Second, we assume 
that the element abundance due to nuclear reaction does 
not differ significantly from the initial abundance. This is 
obviously not a robust assumption if significant nucleosyn- 
thesis does occur in the accretion flow. Before we know 
the final result, we must however make this assumption. 
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We could complete some "iteration" if we knew how the 
abundance changes. Third, in the large-scale universe, the 
relative abundance of an isotope compared to other iso- 
topes of the same element is quite constant. You can find 
the relative abundance of isotopes from the table at the 
website 0. In our calculation here, we accept the assume 
isotope abundances given in this table. 

The nuclear reactions in the star and in the hot accre- 
tion flow are different. In the star, different sets of burn- 
ings are well separated by appreciable temperature dif- 
ferences in the star. However, in the hot accretion flow, 
these nuclear reactions can be ignited almost simultane- 
ously because the temperature is sufficiently high to ignite 
different nuclear burning at the same time. In principle, 
nuclear reaction can proceed for many kinds of elements if 
the temperature is high enough. We considered only H, He, 
and C burning after finding that other nuclear reactions 
were not important contributors to the energy generation 
rate. Our second assumption is therefore justified. 

The temperature of the accretion flow increases when 
it is closer to the central black hole, when an increasing 
number of elements are burned. Different nuclear reactions 
are ignited at different temperatures. Hydrogen is burned 
when the temperature is higher than 7 x 10 6 K. Helium 
can be burned when the temperature reaches 1 — 2 x 10 8 
K, and for carbon, this temperature is about 7 x 10 8 K. 
The total energy generation rates in different temperature 
ranges are as follows: 



(He, 

(He + £CS 



T > 7 x 10 6 K, 
T > 2 x 10 a K, 
T > 7 x 10 8 K, 



(6) 



(h j (He , and ec is the total energy generation rate of hy- 
drogen, helium, and carbon, respectively, and e nuc repre- 
sents their sum. They are defined to be the energy pro- 
duced per unit mass and time 
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where ry is the reaction rate per unit i.e., volume and 
time between the nuclei i and j, and Qij is the energy 
released per reaction, 
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(av)ij is the average cross-section, and and rij is the 



number density of particles 
written as 

rii = pN A ^- = pN A Yi, 

Ai 



and j, which can also be 



(10) 



where X^, Y^, and are the mass abundance, number 
abundance, and atomic weight of element i, respectively, 
and N A is Avogadro number. According to Eqs. (7), (8), 
(9) and (10), the energy generation rate of nuclear reac- 
tions can be written as 

pNaY^t^Y^Q^Na^v)^. (11) 



1 + Si- 



The detailed calculations can be found in Huang (1998), 
Fowler et al. (1975), Harris et al. (1983), and Georgeanne 
& Fowler (1988). 

We study the energy generation rates of nuclear re- 
actions for accretion flows of various parameters. For the 
mass of the black hole, we consider a stellar mass black 
hole (IOMq), an intermediate-mass black hole (1O 3 M0), 
and a supermassive black hole (10 8 M Q ). These values are 
typical of black hole X-ray binaries, some possible ultra- 
luminous X-ray sources, and AGNs, respectively. For the 
case of supermassive black holes, we consider various mass 
accretion rates, namely 10~ i M E (Me is the Eddington 
accretion rate, M E = 10L Edd /c 2 ), 10 _3 M B , O.OIMe, and 
0.1Mb- 

After deriving the global solutions of the hot accre- 
tion flow by solving Eqs. (1-5), we measure the physical 
quantities of the hot accretion flow as a function of radius 
in each case, including temperature Tj and T e and mass 
density p, and calculate the viscosity heating rate q + (de- 
fined as the viscosity heating rate per unit volume). We 
then calculate the energy generation rate e„ uc of the nu- 
clear reactions at different radii based on temperatures 
and densities obtained above from Eq. (11). We note that 
tnuc is the energy generation rate per unit mass, but q + 
is usually defined as the viscosity heating rate per unit 
volume. We therefore define £ nuc = pe nuc for consistency, 
and also denote q + by £„j s . Finally, we derive the ratio 
of energy generation rate of nuclear reaction to viscosity 
heating rate, i.e., e nuc /e V i S . If this ratio is high, the ADAF 
model must obviously be revised by including the heating 
caused by nuclear reactions. 

2.3. Calculation results 

The calculation results are shown in Figs. 1 & 2. Figure 
1 shows the ratio of the heating rate for nuclear reactions 
to that for viscous dissipation for M — O.IMe but dif- 
ferent black hole masses. We find from this figure that 
the ratio is roughly proportional to (r/r s )~ 3 . This is be- 
cause in our case s nuc is roughly proportional to p 2 T 5 , 
while e V i S oc pr~ h l 2 ~ s . So e nuc /e V i S oc pr 5 ^ 2+s T 5 . Since 
p oc (r/r s ) -3 / 2+s while T oc r _1 , we have e nU c/£vis oc 
r 2s ~ 4 oc r~ 3 ' 4 . The ratio does not depend on the mass of 
the black hole. 

Figure 2 shows this ratio as a function of radius for 
different M but the same black hole mass. The ratio is 
roughly proportional to M (or density), because the nu- 
clear reaction heating rate is proportional to M 2 while the 



1 http://www.chem.queensu.ca/FACILITIES/NMR/nmr/massviscous heating is proportional to M, and the ions temper- 
spec/mstable3.htm ature is always virial, independent of M. However, even 



4 



Hui Zhang et al.: Is the energy generation rate of nuclear reactions in hot accretion flows important? 



0.01 p 1 — 




le-07 1 1 1 — 1 1 — 1 1 — 

1 10 100 1000 

r/r 

Fig. 1. The ratio of nuclear-reaction energy generation 
rate to viscous heating rate as a function of radius for 
different black hole masses. The solid, dashed, and dot- 
ted lines represent the cases of M = 1OM , 10 3 M Q , and 
10 8 Mq. The mass accretion rate M — 0.1M E in all three 
cases. 
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Fig. 2. The ratio of energy generation rate of nuclear reac- 
tion to the viscous heating rate as a function of radius for 
different mass accretion rates. The solid, dotted, dashed, 
and dot-dashed lines represent the cases of M = O.IMe, 
10- 2 M E , 10~ 3 M E , and 10" 4 M B . The mass of the black 
hole is 10 8 M CT . 



at M — O.IMe, almost the highest accretion rate a hot 
accretion flow can have, this ratio is less than 1%, i.e., the 
heating rate of nuclear reactions in the hot accretion flow 
is not important. We therefore do not need to include the 
nuclear reaction in calculating the dynamics of ADAF. 



In our calculation, the value of the viscosity parameter 
a is adopted to be 0.3, which is "typical" of ADAFs. We 
should emphasize that this value has some uncertainties. 
The value of a is perhaps constrained most reliably by 
observations of dwarf nova outbursts, where a ~ 0.1 — 0.2 
(Smak 1984). Three dimensional magnetohydrodynamic 
(MHD) numerical simulation suggest that a <~ 10~ 2 — 10 _1 
(Hawley k Krolik 2001). King et al. (2007) argued that 
according to the highest quality observational evidence, 
a ^ 0.1 — 0.4 in thin and fully ionized discs. For a given 
accretion rate, the density of the flow is proportional to 
a -1 . The ratio of the heating rate between nuclear reac- 
tion and viscous dissipation is therefore proportional to 



3. Summary and discussion 

The temperature of the hot accretion flow such as ADAFs 
is very high, so nuclear reaction should occur. To study the 
nucleosynthesis, the first step is to determine the dynam- 
ics of the accretion flow. For this purpose, we have investi- 
gated the ratio of the heating rate for the nuclear reactions 
to that for the viscous dissipation. We have found that the 
ratio increases with accretion rates. However, even for the 
highest accretion rate for which the hot accretion flow so- 
lution exists, the ratio is <, 1%. Therefore, the heating 
caused by nuclear reaction can be neglected. This result 
indicates that the previously obtained global solution of 
hot accretion flow is self-consistent and can be used di- 
rectly in future studies of nucleosynthesis. 

As we mentioned in introduction, some authors 
have previously studied nuclear reactions and nucle- 
osynthesis in accretion flows around black holes (e.g., 
Jin et al. 1989, Arai & Hashimoto 1992, Chakrabarti 
& Mukhopadhyay 1999, Mukhopadhyay & Chakrabarti 
2000, Mukhopadhyay & Chakrabarti 2001). They found 
that significant nuclear reactions can occur and many 
heavy elements can be synthesized, which is quite dif- 
ferent from our results. There are three reasons for this 
discrepancy The first is that, as we mentioned in the in- 
troduction, the works by Jin et al. and Arai & Hashimoto 
are based on a "thick disk" model in which the radial ve- 
locity is several orders of magnitude lower than that in 
ADAFs if the accretion rate is the same. Secondly, the 
accretion rates in most of the above-mentioned works are 
super-Eddington, while the accretion rate in ADAFs is 
lower than the Eddington rate. The third reason is that 
the value of the viscous parameter adopted in previous 
work is much smaller than that we adopt, a = 10~ 7 
(Jin et al. 1989), 10~ 10 (Arai & Hashimoto 1992), and 
(Chakrabarti & Mukhopadhyay 1999; Mukhopadhyay & 
Chakrabarti 2000; Mukhopadhyay & Chakrabarti 2001), 
respectively. So the density of the accretion flow will be 
much higher than that in ADAFs. Since the efficiency of 
nuclear synthesis will increase in proportion to density, 
this explains why they found that the nuclear reaction is 
so significant. 
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